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The  Angular  Impact  Distibution  of 
Charged  Particles  Attracted  to  a 
Charged  Cylindrical  Spacecraft 


1.  INTRODUCTION 

Spacecraft  in  geosynchronous  orbit  exhibit  charging  when  immersed  in  hot 
substorm  plasmas.  The  potential  that  a  surf-ice  attains  depends  both  on  the 
material's  properties  and  on  the  plasma's  characteristics.  Secondary  emission 
of  electrons  by  electrons  and  ions,  as  well  as  backsrattering,  determine  the 
floating  potential.  These  are  both  functions  of  the  angle  of  incidence  of  particles 
on  the  surface. 

In  this  work,  the  angular  distribution  of  incident  particles  is  obtained  for 
attracted  and  repelled  particles  for  cylindrical  spacecraft.  Results  are  given 
for  the  plane  and  spherical  cases  as  well. 

For  isotropic  spacecraft  materials,  the  backscattering  coefficients,  the 
secondary  emission  coefficients,  and  the  shielding  efficiency  depend  on  impact 
energy  and  impact  zenith  angle,  but  are  independent  of  azimuth  angle.  According¬ 
ly,  only  the  distribution  of  impact  zenith  angles  is  important.  For  convenience, 
only  protons  will  be  considered,  the  extension  to  other  charged  particles  being 
obvious.  First  the  problem  will  be  precisely  defined,  then  a  simplified  model  of 
physical  reality  set  up,  and,  finally,  the  distribution  of  impact  angles  derived. 
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2.  DK.KIMTION  Oh  PROBLEM 


Kind  the  flux,  j  of  protons  striking  an  element  of  area  AA  with  impact 
kinetic  energies  within  the  band  K  t  Jlv'2  and  with  impact  zenith  tingles  within 
the  band  0  t  A??/2. 


3.  MODEL 


The  model  is  that  of  an  infinitivelv  long  cylinder  of  radius  X  .  Revotid  a 

l  ) 

"sheath"  of  radius  R^,  the  plasma  is  isotropic  and  is  at  zero  potential.  The 
potential  within  the  sheath  is  modeled  by: 


ln(K/R  .) 

o  ln(R  /R  ) 
o  s 


(1) 


1'he  subscripts  "o"  and  "s"  will  be  used  to  indicate  values  of  variables  at  R  and 

R  ..  No  subscript  will  be  used  for  general  values  of  R. 

&  -■>  -1-1-1 
In  the  model,  the  differential  flux  f(K,  protons  cm  “  see  ster  eV  is 

a  function  of  kinetic  and  potential  energies  oniy  —  except  in  the  "shadow"  of  the 

Satellite  where  it  is  zero.  Both  potential  ($)  and  kinetic  (K)  energy  are  expressed 

in  electron  volts.  Any  impact  angle  and  energy  corresponding  to  a  trajectory 

originating  elsewhere  on  the  satellite  is  considered  to  be  "shadowed.” 


4.  DERIVATION  OK  IMPACT  ZENITH  ANGLE  DISTRIBUTION 

The  impact  proton  flux  is  the  product  of  four  factors,  namely:  energy  band 
width,  unshadowed  solid  angle,  projected  area,  and  differential  flux.  Tins  may 
be  w  ritten  as 

j  -  (AEKA0  d>  sin  0  )<AA  cos  0  )[f(E  ,  «  ))  (2) 

•O  7.  O  O  O’  O 

where  =  total  unshadowed  azimuth  angle.  The  solid  angle  is  that  which  would 

be  subtended  at  the  earth's  center  by  the  area  enclosed  by  latitudes  D  ±  At?  2  and 

‘  o 

longitudes  ij>  ±  $z/2.  The  azimuth  angle  will  be  taken  as  zero  perpendicular  to  the 
cylinder's  axis.  If  shadowing  occurs,  it  will  be  at  the  larger  angles. 

Equation  (2)  may  be  rewritten  as 

j  ;  (AE  A<?  AA  f(E  ,  9  )1  [sin  8  cos  B  )  (4*  )  (3) 

Jo  o  o  o  o  ' 


fi 


where  bm  -  maximum  unshadowed  azimuth  angle,  <t>  being  measured  So  as  always 
to  lie  in  the  first  quadrant.  The  factor  4  accounts  for  the  four  quadrants.  A 
further  rewriting  yields 

i  =  K  <f>  sin  2  0 
■o  tn  o 

K  =  ^  AO  AA  f(E  ,  $  )  (4) 

2  o  o 

0  <  0  c  ff  /  2  . 

““  m 

Note:  sin  2  0  =2  sin  0  c os  0 
o  o  o 

Equation  (4)  is  very  general.  The  geometry  of  the  spacecraft's  surface  and  field 
enter  into  the  equation  only  through  $  the  maximum  unshadowed  azimuth  angle. 

5.  SPECIAL  CASES 

For  positive  spacecraft  potentials,  the  reverse  trajectories  curve  away  from 
the  spacecraft;  therefore,  there  is  no  shadowing  and 

<1>  ---  v/2  if  4  ?  0  .  (a) 

From  Eq.  (4),  the  angular  distribution  is; 

j  =  (tr/2)  K  sin  2  0  . 

■  o  ' 

The  other  special  case  is  that  of  looking  for  protons  of  impossibly  low  energy. 
Here  we  may  write 

<j>  -  0  if  9  <  -E  ,  ((>) 

m  o  o 

which  has  meaning  only  for  negatively  charged  spacecraft. 


6.  MAXIMUM  UNSHADOWED  AZIMUTH  ANGLE 

The  simplest  cases  have  been  dealt  with  in  the  preceding  section.  This 
section  will  consider  cases  defined  by: 


(7) 


-E  <  $>  <0 

o  o 


The  model  has  zero  axial  field.  Therefore,  the  kinetic  energy  associated  with 
the  axial  velocity  component  is  independent  of  radial  position  and  may  be  written: 

E  -  E  sin2  <j>  sin2  6  ,  (8) 

a  o  o  o 

where  E  is  the  impact  energy  and  <4o  and  0^  are  the  impact  angles.  This  rela¬ 
tionship  follows  from  the  breakdown  of  the  impact  velocity  into  four  components, 
as  shown  in  Figure  1.  The  model  has  a  radial  field  about  the  cylinder's  axis. 
Thus  the  tangential  (l  to  radial  and  axial  directions)  kinetic  energy  is  subject  to 
conservation  of  angular  momentum  and  is  given  by: 

E  (R)  =  E  (R  /R)2  cos2  <4  sin2  6  .  (9) 

t  o  o  o  o 

The  radial  kinetic  energy  is  simply: 


Ep  =  E  cos"  0  .  (10) 

The  total  proton  energy  is  independent  of  the  radius,  and  beyond  the  sheath  it 
equals  kinetic  energy.  The  total  energy  is  thus: 

E  =  E+  E+Ef  +  $  ,  <  0  .  (11) 

rat 

At  the  sheath's  edge  Eq.  (11)  becomes: 

E  =  E  cos2  0  +  E  sin2  <4  sin2  0  +  (R  /R  )2  E  cos2  <4  sin2  0  (12) 

s  so  o  o  o'  S  o  o  o 

and  at  the  satellite's  surface  it  becomes 

E  =  E  cos2  0  +  E  sin2  $  sin2  0  +  E  cos2  <4  sin2  0  +  $  (13a! 

o  oo  o  oo  o  oo 


E  =  E  +  $ 
s  o 


$  <  0  . 


The  limiting  case  corresponds  to  a  vanishing  radial  energy  at  the  sheath's  edge  or 


cos  0=0 
s 


I 
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Equations  (4)  and  (lt>)  together  give  the  zenith  angle  (0)  distribution  of  protons 
striking  the  surface  of  a  cylindrical  satellite. 

The  angular  distribution  derivation  given  has  considered  only  the  potential  at 
the  inner  (R  -  R  )  and  outer  (R  =  Rg)  edges,  not  the  potential  within  the  sheath 
(R  <  R  ,).  It  may  be  shown  that  the  potential  within  the  sheath  is  not  controlling, 

s  '  i 

provided  that  at  all  points  within  the  sheath  it  is  equal  to  or  more  negative  than 
that  of  the  inverse  square  potential  function: 


$ 


(R  /R)2- 1 
s _ 

(r ,/r  r- i 


<t  <  0 

o 


(17) 


This  function  is  adjusted  to  match  the  model  at  the  sheath's  edge.  The  logarith¬ 
mic  function  of  the  model  satisfies  this  condition.  The  proof  is  similar  to  that 
given  for  the  spherical  case  for  objects  in  a  gravitational  well.  The  proof  is  only 
valid  if  the  potential  is  a  function  of  radial  position  alone. 

The  distribution  of  impact  zenith  angles  according  to  Eqs.  (4)  and  (lii)  is 
shown  in  Figure  2  for  several  cases.  The  deviations  from  the  symmetrical 
sin  (2  0)  distribution  are  due  entirely  to  shadowing.  The  shadowing  is  equivalent 
to  a  "horizon"  at  less  than  90°  below  the  zenith  due  to  refraction  of  the  protons. 


0°  10°  20°  30°  40°  50°  60°  70°  80°  90° 


ZENITH  ANGLE 

Figure  2.  Distribution  of  Zenith  Impact  Angles  for  Protons  on  a  Cylindrical 
Satellite:  (a)  Thin-Sheath  C  ase 
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T.  DISCI  SS|()N 


I'hi'  •hf.iy  given  here  is  exact  in  i  licit  it  is  i  .  •  mph  ely  n  msis'  •  oh  the 
model.  The  model,  hi nvrver,  r epresents  a  physical  impossibility  in  that  t li«  *-!•*»  - 
•  ric  field  is  discontinuous  tit  flic  surface  of  the  sheath  and  an  approximation  ;n 
that  end  effects  are  taken  to  be  vanishingly  smalt.  Only  if  the  nearest  end  is  tit  a 
distance  very  large  compared  to  the  radius  will  its  effect  become  vanishingly 
small.  file  theory  does,  however,  represent  a  limiting  case  and,  as  such,  is 
useful. 


8.  SIMM  VIO 


The  angular  distribution  of  impacting  repelled  particles  is  sin  (2  '>),  as  given 
by  Kq.  (4),  since  there  is  no  shadowing  for  repelled  particles,  electrons  inci¬ 
dent  on  a  negatively  charged  surface  and  ions  on  a  positively  charged  surface  w  ill 
impact  in  this  way.  In  order  to  calculate  surface  potentials,  taking  account  "f 
sec  ondary  electron  emission,  one  must  know  the  distribution  of  impact  angles  a., 
well  as  the  secondary  electron  yield  as  a  function  of  tingle. 

The  angular  distribution  of  impacting  attracted  particles  (electrons  <>n  a 
posif’ve  surface,  and  positive  ions  on  a  negative  surface)  is  more  complicated. 

The  impacting  angular  distribution  fin-  a  cylindrical  surface  for  attracted  parti.  1«  s 
has  been  derived  here,  as  shown  in  Figures  2  (a),  2  (b),  and  2  (<•),  for  a  range-  of 
energies  and  sheath  thickness. 

The  distribution  of  impact  angles  (attracted  particles)  for  a  plane  surface  is 

Mi*  same  as  'he  case  of  zero  sheath  thickness.  For  this  rase  the  distribution  is 

proportional  to  sin  (2  0)  for  0  less  than  0  ,  where 

m 

-  cos'1  (4>/K)1/2  , 

and  where  K  is  the  impact  energy.  If  6  is  greater  than  0  the  distribution  func¬ 
tion  is  zero.  For  spherical  surfaces  there  is  no  shadowing,  so  the  distribution  is 
again  proportional  to  sin  2  0  in  the  thick  sheath  limit. 
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